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Abstract 

The production of charged pions in minimum bias p+C interactions is studied using a 
sample of 377 000 inelastic events obtained with the NA49 detector at the CERN SPS 
at 158 GeV/c beam momentum. The data cover a phase space area ranging from to 
1.8 GeV/c in transverse momentum and from -0.1 to 0.5 in Feynman x. Inclusive invariant 
cross sections are given on a grid of 270 bins per charge thus offering for the first time 
a dense coverage of the projectile hemisphere and of the cross-over region into the target 
fragmentation zone. 



*^ Corresponding author: Siegfried. Wenig@cern.cli 



1 Introduction 



The study of p+A interactions represents an important part of the NA49 experimental 
programme which is aimed at a comprehensive study of soft hadronic interactions covering 
both elementary and nuclear collisions. In the framework of this programme a first paper on 
inclusive pion production in p+p collisions has been published recently iUl. The present study 
extends these measurements to p+C interactions. The use of the isoscalar, light ^^C nucleus 
has several motivations. Firstly, it allows inspection of the evolution from elementary to nu- 
clear collisions for a small number of intranuclear projectile interactions as compared to the 
high-statistics data sets on the heavy ^o'^Pb nucleus from the NA49 experiment [2J. Secondly, it 
satisfies, at least partially, the need of high precision p-i-A reference data from light nuclei for 
the control of systematic effects in neutrino physics. This applies both to long baseline neutrino 
oscillation experiments using light neutrino production targets (e.g. |3 1 which uses a beam en- 
ergy of 120 GeV) and to atmospheric neutrino studies [4J. In the latter case the ^^C nucleus is 
sufficiently close to ^^N and ^^O to allow a sensible test of the different hadronic production 
models developed for the interpretation of these data [5|. 

The NA49 p+C data sample consists of 377 000 inelastic events corresponding to only 
8% of the event number available in p-i-p collisions 1 1 1. The reasons for this relatively limited 
statistics lie entirely in the heavy restrictions which the SPS fixed-target programme as a whole 
has suffered over the past years. These restrictions concern both the beam availability and the 
choice of topics in hadronic physics acceptable to the relevant committee |l6|. In view of the al- 
most total absence of existing data in the SPS energy range CHHI the NA49 results nevertheless 
constitute a unique extension of the present knowledge in the sector of p-i-A interactions with 
light nuclei. 

The present study is closely related to the preceding publication on pion production in 
p+p interactions [ 1 1 which is recommended as a reference for the description of most exper- 
imental details. Only those items which are specific for p+C collisions will be touched upon 
here. The layout of the paper is arranged as follows. After a short comment concerning the sit- 
uation of previous measured data in Sect. 2, Sect. 3 describes the parts of the NA49 experiment 
which differ from the p+p data taking. Sect. 4 deals with the determination of the inclusive 
cross sections and the applied corrections. The final data are presented and compared to other 
experiments in Sects. 5 and 6. The pt integrated distributions are evaluated for minimum bias 
condition in Sect. 7 and Sect. 8 describes the dependence on the measured number of grey 
protons. 

2 The Experimental Situation 

The asymmetric nature of p+A interactions necessitates in principle a complete coverage 
of the target and projectile hemispheres in order to experimentally constrain the underlying pro- 
duction mechanism. The phase space coverage of exisiting data in p+C is even more restricted 
than for p+p collisions [1|. In fact there are only two sets of data in the SPS energy range 
offering double differential inclusive cross sections for identified pions. 



where the longitudinal scaling variable xp = 2pl/ ^ is defined in the nucleon-nucleon cms. 

The first data set at 400 GeV/c beam momentum covers the far backward region at a 
number of fixed laboratory angles [8|, yielding important information about the region of intra- 
nuclear cascading. The second one presents a small number of data points in the far forward 
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Figure 1 : Phase space coverage of existing data in a) Fermilab and b) NA49 



region Q. The situation is depicted in Fig.[Tt which shows that there are no data in the complete 
central region -0.2<xf<0.3. 

The new NA49 results fill this gap to a large extent, as shown in Fig.[Tj5. They are limited 
to < 1.8 GeV/c and to xp < 0.5 essentially by the modest event statistics, see Sect.[ll In the 
backward hemisphere, a px dependent cut atxp > -0. 1 is imposed both by the NA49 acceptance 
and by the particle identification via ionization energy loss. Nevertheless an inspection of the 
important cross-over region between projectile and target hemispheres is possible. 



3 Parts of the NA49 experiment specific to p+C interactions 

Most parts of the NA49 experimental setup and of the data processing procedure are 
identical for p-i-p and p-i-A data taking. This concerns the beam definition, the interaction trig- 
ger scheme, tracking, event reconstruction and selection, and particle identification [1|. The 
experimental items which are specific to the p-i-C data taking are described below. 

3.1 Target and grey proton detection 

A graphite target of 0.7 cm length and 0.6 cm diameter with a density of 1.83 g/cm^ has 
been used. This corresponds to an interaction length of 1 .5%. The target is housed inside a grey 
proton detector (9 1 which has been developed for the control of impact parameter (centrality) in 
p-i-A collisions and is presented schematically in Fig.|2l 

This detector is a cylindrical proportional counter of 12 cm diameter which surrounds 
the target and has a window in the forward hemisphere corresponding to the acceptance of 
the spectrometer inside polar angles of < 45°. Its surface is subdivided into 256 pads which 
provides ample granularity for the counting of the typically less than 8 grey protons measured 
per event in light ion applications. A thin (200 /xm) copper sheet on the inner surface absorbs 
nuclear fragments by range, and an electronics threshold placed at 1.5 times the minimum ion- 
ization deposit cuts high momentum particles as the grey protons are placed high on the 
branch of the Bethe-Bloch energy-loss distribution. Grey protons in the momentum range 0.15 
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Figure 2: Centrality detector 



to 1.2 GeV/c reconstructed and identified inside the spectrometer acceptance are added to the 
number measured in the centrality detector. 

The resulting number distribution is shown in Fig. |3l together with the prediction by 
Hegab and Hiifner [ 10 1 for the total grey proton yield and the measurement by Braune et al. 
fl 1 1 at 50 and 100 GeV/c beam momentum at the SPS using a detector of 60% geometrical 
acceptance. Compared to this the NA49 data correspond to a 45% effective acceptance which 
is in part due to the steep angular distribution of slow protons in p+C collisions ifTDl and to the 
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Figure 3: Number of grey protons rigj-ey distribution. The lines through the data points are drawn 
to guide the eye 
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large opening of ±45° imposed by the spectrometer acceptance. Only about 30% of protons in 
this angular range with momenta below 1 GeV/c are in fact reconstructed and identified by the 
NA49 tracking system. 

Due to the short available data taking period and to the sharp drop of the event yield as 
a function of the number of grey protons (less than 30% of all events have a measured grey 
proton), the on-line triggering capability on grey protons available in the NA49 trigger system 
could not be used. All data have therefore been obtained in "minimum bias" condition without 
imposing centrality selection. It is nevertheless possible to use the grey proton information in a 
sample of px integrated inclusive cross sections in order to study the evolution with centrality 
as described in Sect.lU 

As the range of grey protons at their most probable momentum of 0.3 GeV/c is only 
1 cm in Carbon, the target diameter of 6 mm has been kept at a minimum with respect to the 
transverse beam profile [ 1 1 in order to reduce the absorption of grey protons by energy loss 
in the target material. The small fraction of beam particles in the tail of the transverse profile 
beyond the target radius 3%) has been corrected for in the determination of the trigger cross 
section. 



3.2 Trigger cross section 

The trigger scheme was the same as in p-i-p interactions, using an interaction trigger de- 
fined by a small scintillation counter 380 cm downstream of the target in anti-coincidence with 
the beam. Due to the reduction of forward protons in p+A interactions by baryon number trans- 
fer towards central rapidities and due to the corresponding yield decrease of produced particles 
at large xp, the systematic effects connected to this trigger method are smaller than in p-i-p colli- 
sions and amount to a reduction of the trigger cross section of 9%±2% with respect to the total 
inelastic cross section. A break down of these losses determined by using measured inclusive 
distributions of protons, pions and kaons [ 17 1 in a Monte Carlo calculation, is given in Table[T] 



0"trig 

loss from p 
loss from tt and K 
contribution from cTgi 
predicted a^nci 


210.1 mb 
17.1 mb 
2.4 mb 
-3.3 mb 

226.3 mb 


literature value 


225.8 mb 



Table 1 : Contributions derived from the detailed Monte Carlo calculation to the determination 
of the measured inelastic cross section cTinei 

The resulting measured value of the inelastic p+C cross section is 226.3 mb with an esti- 
mated systematic error of 2.5% (see Table|3l). It compares well with a compilation of preceding 
measurements ifOI giving 225.8 mb as shown in Fig. HI 

The relatively large systematic deviations visible in the reference values, which have all 
been obtained in transmission experiments using yield extrapolation to zero momentum transfer, 
is noteworthy. It is partially attributable to the definition of the elastic component and to the 
corresponding uncertainty in the slopes of the t-distributions used. Compared to the precision 
available in the elementary hadron-nucleon inelastic cross sections, the absolute normalization 
of inclusive yields therefore suffers from larger systematic uncertainties. 
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Figure 4: Inelastic cross section compared to previous measurements ifTSi 

3.3 Event sample and cuts 

In addition to the cuts imposed on the transverse beam definition as obtained from the 
Beam Position Detectors [ 9 1 an additional cut on the transverse beam radius at less than 3 mm is 
imposed due to the target dimension (Sect. lTTT) . The longitudinal vertex position is constrained 
to a fiducial region, as exemplified for two values of charged multiplicity in Fig.|5l where the 
distribution of detector material in the vicinity of the target is clearly visible. 

The combined cuts result in the event sample given in Table |2l where in particular the 





Figure 5: Normalized vertex distributions from full and empty target events with selected track 
multiplicity a) one and b) five and more. The arrows indicate the fiducial region 
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reduction of the empty target fraction from 30% in the total sample to 16% after cuts is note- 
worthy. Due to the small correction imposed by the empty target contribution (see Sect. 14. IL 
the fraction of running time spent on empty target could be kept at about 5% only. 



Events taken 
Full target Empty target 


Events after selection 
Full target Empty target 


535.7 k 31.2 k 


377.6 k 11.8 k 



Table 2: Data sample 



3.4 Acceptance coverage, binning and statistical errors 

Due to the forward-backward asymmetry of p-i-A interactions a substantial effort was 
spent in exploiting the available acceptance in the region of negative xp- Here limitation are 
imposed both by the detector acceptance itself and by the particle identification problems in the 
region near minimum ionization. The resulting phase space coverage with 270 bins per charge 
is presented in Fig. |6l where also the statistical errors per bin are indicated. 

The reduced statistics compared to the much larger event sample in p-i-p collisions Q 
results in larger bins generally and a limitation to the regions < 1.8 GeV/c and xp < 0.5. It 
nevertheless provides an unprecedented overall coverage which allows a detailed study of the 
important cross-over region between target and projectile fragmentation for the first time. 




Figure 6: Binning scheme inxp and pt together with information of the statistical error 



3.5 Particle identification 



The procedures for the extraction of pion yields from the energy loss distributions mea- 
sured in each bin are identical to the ones used in p+p collisions [IJ for most of the forward 
hemisphere. The extended coverage of the backward hemisphere together with the forward- 
backward asymmetry of p-i-C interactions necessitates an extension of the methods developed 
in im for the treatment of the lab momentum region below 3 GeV/c where the energy loss func- 
tions of pions, kaons and protons approach each other. The kinematic situation in this region is 
indicated in Fig.|71 




Figure 7: dE/dx cross-over region: a) in {xf,Pt) plane of the pions, b) in {xp,Xp) plane, and c) 
in (xp,Xp) plane. The unhatched areas correspond to the regions where pions can be identified 
due to the increase of energy loss of the kaons and protons, respectively 

Using the pion mass in the transformation from xp to lab momentum the critical zone is 
defined by the hatched area in Fig.|7^ where the upper limit corresponds to pi^h = 3 GeV/c and 
the lower limit traces the px cut-off used in the data extraction (Fig. |6l). This area is mapped 
into the xp regions for kaons (Fig.lTj?) and protons (Fig.|7|:) as a function of a;^°" when using 
the proper masses in the corresponding Lorentz transformations. 

For protons the lower part of the critical zone is again available for dE/dx extraction due 
to the rapid increase of their energy loss in the region of the Bethe-Bloch function. The 
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remaining band at -0.4 < < -0.2 is treated by reflecting the xp bins into the forward 

hemisphere and by using an interpolation of the measured backward/forward yield ratios shown 
in Fig. H 

The p/ 7r~ and K/tt ratios are smaller than 10% for all bins in the critical area. Here again 
the method of bin reflection was used by imposing an extrapolation of the measured yield ratio 
with respect to p+p interactions into the backward hemisphere as shown in Fig. 15] for kaons. 
In this case the ratio in the far backward region corresponds to the mean number of projectile 
collisions predicted from the inelastic p-i-C and p+p cross sections ifT^ . For kaons the different 
loss rates from weak decays in the forward and backward hemispheres were taken into account 
in the yield determination. 

4 Evaluation of invariant cross sections and corrections 

The invariant inclusive cross section 

(Pa 

f{xF,PT) = E{xf,Pt) ■ -r^{xF,PT) (2) 
is experimentally defined by the measured quantity Q 

. / A 3\ Tpf A 3N ^trig An{xF,PT,Ap^) 

/meas(a;F,PT, ^P ) = E{xf,Pt, ^P ) " ^TT 7"^ , (3) 

where Ap^ is the finite phase space element defined by the bin width. 

As in lUl several steps of normalization and correction have to be applied in order to make 
fmeasi^F^PT, ^P^) approach f(xF,PT)- The determination of the trigger cross section and its 
deviation from the total inelastic cross section has been discussed above. The corrections for 
pion weak decay and absorption in the detector material are identical to p+p and will not be 
discussed here. The remaining corrections which are numerically different in p+C interactions 
are: 

- treatment of empty target contribution 

- re-interaction in the target volume 

- effect of the interaction trigger 
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- feed-down from weak decays of strange particles 

- effect of finite bin width. 

These corrections will be described and quantified below. 

4.1 Empty target contribution 

In the determination of the normalized quantity 




the ratio e of empty over full trigger rates has been reduced from 0.3 to 0. 1 6 by the cuts described 
in Sect. l3.3l above. In addition, the bulk of the remaining empty target yield is produced in Mylar 
foils and air (Fig.|5l). These materials are sufficiently close to Carbon to make the normalized 
bin contents (An/Ne^)^ and (An/iVev)^^ approximately equal. The deviation of the complete 
normalized yield from the full target yield alone is therefore expected to be small and to be 
essentially defined by the different fraction of empty events in full and empty target conditions. 
This is demonstrated in Fig. flOl where the ratio between the bin contents for full-empty target 
and full target alone is shown as a function of xp. 
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Figure 10: Correction factor applied to account for empty target contribution as a function of 
xp for the average of 7r+ and tt^ 

Within the limits of the experimental accuracy no differences between 7r+ and vr and no 
dependencies on px have been observed. 



4.2 Target re-interaction 

The Carbon target has an interaction length of 1 .5% which corresponds to only 55% of the 
length of the hydrogen target used in p-i-p interactions. The expected re-interaction correction 
is therefore below the 2% level even in the extreme backwards bins and has been scaled down 
accordingly from the values obtained in H). 
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4.3 Trigger bias correction 



Several effects contribute to a modification of the correction for the trigger bias which is 
introduced by the interaction trigger in p+C as compared to p+p collisions. They all lead to a 
reduction of the correction to the p+C data. 

- Due to baryon number transfer towards central rapidity ("stopping") there are less for- 
ward protons hence a smaller probability to veto events by the trigger counter. 

- Due to the correlated steepening of the x p distributions of produced particles (see Sect.|5l) 
there is again a reduction of the veto probability. These two effects combine to explain the 
higher inelastic trigger efficiency of 93% in p+C as compared to 89% in p+p collisions. 

- Unlike in p+p events, there is a long-range correlation between target fragmentation and 
forward particle density. Large xp protons are correlated with single projectile collisions 
yielding small backward multiplicities, whereas multiple projectile collisions result in 
high target multiplicity and low forward yields. This correlation reduces the effect of the 
interaction trigger in the backward hemisphere. 

The trigger bias correction has been obtained using the method developed in [ 1 1 by artifi- 
cially increasing the diameter of the trigger counter in the analysis and extrapolating to surface 
zero. The resulting corrections are shown in Fig.[TT]as a function of for two values of px- 
The correction is smaller than in the p+p case lU and thus confirms the reduction quoted above. 
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Figure 11: Trigger bias correction as a function of xj? at various pp for a) 7r+ and b) tt" 



4.4 Feed-down correction 

A principle problem in the determination of the feed-down correction induced by the 
weak decays of strange particles lies in the absence of data on and strange baryon production 
in light ion collisions. The corresponding yields have therefore been determined from the NA49 
data directly using the following yield ratios with respect to p+p interactions: 

- The yield is extracted from (K+ + K^)/2. The measured ratio is extrapolated into 
the backward hemisphere using a two-component superposition picture L2J as shown in 
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Figure 12: Yield ratios for a) K^, b) A and c) A 



Fig. fT2k . As the kaon yields do not suffer from isospin effects f?] this extrapolation is 
straight-forward concerning the target contribution. 
- The evolution of A and A yields relative to p+p is obtained from p+vr" and p+7r+ mass dis- 
tributions exploiting the event mixing technique described in ifTSi and using vertex tracks 
both for the baryon and for the meson involved. The resulting yield ratios are shown in 
Figs.[T2h and c. For the extrapolation into the backward hemisphere the two-component 
superposition picture is again used. No isospin effects are present in this extrapolation. 

The yield ratios for are derived from the A parametrization using the S/ A ratios from 
p-i-p in the projectile hemisphere. In the target fragmentation region the expected isospin effects 
are taken into account. 
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Figure 13: Feed-down correction to pion from weak decays for a) 7?+ and b) vr 
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The resulting overall feed-down corrections are shown in Fig. [131 as a function of xf for 
several values of pt- 



4.5 Binning correction 



The correction for finite bin width follows the scheme developed in [T| determining the 
deviation of the real cross section at the bin center from the measured one (averaged over the 
bin) using the local second derivative of the particle density distribution. The derivative is cal- 
culated using the experimental data and therefore does not depend on specific parametrizations. 
Typical values of the resulting correction are shown in Fig.[l4|for low px as a function oi xp 
and for xp = ^ function of pt, both for constant bin widths of Axi? = 0.02 and Apj- = 
0. 1 GeV/c, and for the bin widths actually used in the data extraction. 

Due to the increased average bin width in p+C the values are somewhat larger than in 
p+p, but typically stay below the ±2% limit, except for xp > 0.4 and > 1-2 GeV/c. 
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Figure 14: Correction due to the binning in a) xj? and b) p^- The crosses represent the correction 
at fixed bin widths of Axf = 0.02 and Apr = 0.1 GeV/c, respectively and the open circles 
describe the correction for the bins actually used 



4.6 Systematic errors 

An estimation of the systematic errors induced by the overall normalization and by the 
applied corrections is given in Tabled 

An upper limit of 7.5% from linear addition of the error sources and an rms error of 3.8% 
from quadratic summation are obtained. Further information on the bin-by-bin variation of the 
different error contributions is contained in Fig. [15] which demonstrates that the single errors 
fluctuate around well-defined mean values with limited skewness and tails. 

5 Results 

The set of double differential invariant cross sections obtained from the data analysis and 
correction procedures described above forms, by its dense coverage of the available phase space 
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Normalization 


2.5% 


Tracking efficiency 


0.5% 


Trigger bias 


1% 


Feed-down 


1-2.5% 


Detector absorption 




Pion decay vr ^ /i + z/^ 


0.5% 


Re-interaction in the target 




Binning 


0.5% 


Total(upper limit) 


7.5% 


Total(quadratic sum) 


3.8% 



Table 3: Summary of systematic errors 
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Figure 15: Distribution of correction for a) target re-interaction, b) trigger bias, c) absorption in 
detector material, d) pion decay, e) empty target contribution, f) feed-down, g) binning, and h) 
resulting total correction 



in 270 bins per charge, an internally consistent ensemble that reveals, as in the case of p-i-p 
interactions, complex structures which pervade both transverse and longitudinal momentum 
dependencies. These structures defy simple parametrization with straight-forward arithmetic 
expressions. To make full use of the consistency of the data set and to allow for optimum 
precision in the determination of integrated cross sections, a numerical interpolation scheme 
has therefore been used that relies on local continuity in both kinematic variables and allows 
for limited extrapolation into the inaccessible regions of phase space. This chapter summarizes 
the numerical information in data tables, gives a set of distributions as a function of pr, xp and 
rapidity y, and shows the comparison to existing data. 
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5.1 Data tables, distributions and interpolation scheme 

Tables m and 121 present the invariant inclusive cross sections for 7r+ and vr^ respectively. 
They correspond to the binning scheme discussed in Sect. 13 .41 above and reflect the attempt to 
cover the kinematic plane as completely as possible given the limited statistical accuracy of the 
data sample. 

The distributions of the data in pt and xp are shown in Figs.fTTIandfTSlrespectivelv. Here 
the full lines represent the data interpolation mentioned above. 

In order to demonstrate the statistical consistency of the interpolation scheme, the distri- 
bution of the differences between data and interpolation, normalized with the statistical error in 
each data point, is shown in Fig. [161 The shape and width of the distribution comply with the 
expected Gaussian behaviour. 
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Figure 16: Histogram of the differences A between the measured invariant cross sections and 
the corresponding interpolated values (7r+ and tt^ combined) divided by the experimental un- 
certainty A/ of the data points 

5.2 tt^/tt" ratios 

The dependence of the /n^ ratio on and x p carries important information concern- 
ing the details of charge conservation in the hadronization process. As presented in Figs. [T^ and 
EDI important substructures in both variables become visible which are directly comparable to 
the situation in p+p interactions yj. 
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Table 4: Double differential invariant cross section /{xftPt) [mb/(GeV^/c^)] for tt"*" in p+C 
interactions at 158 GeV/c. The statistical uncertainty A/ is given in % 
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Table 5: Double differential invariant cross section /{xftPt) [mb/(GeV^/c^)] for tt in p+C 
interactions at 158 GeV/c. The statistical uncertainty A/ is given in % 
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Figure 17: Invariant cross section as a function of px at fixed xp for a) 7r+ and b) tt produced 
in p+C collisions at 158 GeV/c 
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Figure 19: Ratio of invariant cross section for 7r+ and vr as a function of pr at fixed xp- The 
lines represent the result of the data interpolation 
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Figure 20: Ratio of invariant cross section for 7r+ and vr as a function of xp at fixed px 



This comparison shows global similarities. Remarkable differences are however visible 
in the backward direction where the ratios are expected to approach unity due to the isoscalar 
nature of the target nucleus, and in the detailed structures in the projectile hemisphere. 



5.3 Rapidity and transverse mass distributions 

The rapidity distributions at fixed pt presented in Fig. |2T] extend up to 1.5 units into the 
target hemisphere at px < 0.4 GeV/c and therefore allow a clear view of the asymmetry which 
reaches about 0.25 units at low pt- 

In comparison to p+p collisions [1], they show an important steepening in the projectile 
hemisphere at all transverse momenta, whereas the characteristic deformation at low p^ and y 
is still clearly visible. The mx distributions for both charges are shown in Fig.l22l 

As stressed in 1 1 1 these distributions are all but exponential and the variation of the local 
inverse slopes T with uit, shown in Fig.l^H is even more pronounced than in p+p collisions. 
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Figure 21 : Invariant cross section as a function of y at fixed px for a) 7r+ and b) tt produced in 
p+C collisions at 158 GeV/c 
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Figure 22: Invariant cross section as a function of itit — itIt^ for a) 7r+ and b) tt produced at 
y = 0.0 




Figure 23: Local slope of the mx distribution as a function of — m„ for 7r+ and vr . The line 
shown is to guide the eye 



6 Comparison to other data 

As shown in Sect. |2l the only existing data set which can be directly compared to the 
NA49 results is that of Barton et al.f7|. This comparison is shown in Fig.l24l 

For the 10 points which overlap with the NA49 data a clear upward deviation with an 
average of +25% or 3.6 standard deviations is evident. This deviation is somewhat hard to un- 
derstand as these results come from a group which has published results on p+p interactions 
with the same apparatus lfT4ll which show excellent agreement with NA49 (see lUl for a detailed 
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Figure 24: Comparison of the invariant cross section as a function of xp at fixed px from NA49 
(full symbols) with measurement from |7 1 (open symbols). The n~ data lines are multiplied by 
0. 1 to allow a separation from the 7r+ 

discussion). Also the p+p reference data obtained in the framework of Q are internally consis- 
tent, notwithstanding their sizable statistical errors, both with [ 14 1 and with NA49. Preliminary 
analysis reveals the same problem also with proton yields in p+C collisions. 

As the differences are similar for both pion charges, the /ir^ ratios are expected to be 
unaffected by the problem. This is indeed the case, the good agreement between the two data 
sets at pt = 0.3 GeV/c as a function of Xi? is shown in Fig. |25l 
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Figure 25: Comparison of the tt^/tt ratio as a function of xp at fixed pt from NA49 (full 
circles) with measurement from [j7J (open circles). The dashed line represents the p+p data lUl 

Since the Barton ||7| data extend up to x^? = 0.88 at their measured transverse momenta, 
one can use the consistency of the n^/n^ ratio to extend the NA49 data - albeit with large 
error bars - into the region of large xp- This shows, as indicated in Fig.|25l a sizable increase 
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of the charge ratio with respect to the p+p data also given in the Fig. It is also possible, under 
the assumption that the large systematic yield difference has no dependence on pr or xp, to 
extend the interpolation of the NA49 data towards large xphy imposing a 25% reduction on 
the Barton et al. data. This extension is also indicated in Fig.|24lfor px = 0.3 GeV/c. It shows 
the characteristic break in the xp dependence at between 0.5 and 0.6 also visible in the p+p 
data LU- 

7 Integrated distributions 

The pt integrated yields 



dn/dxp = vr/cTinei • Vs/2 J f / E ■ dp^ 
F= I f-dpl 



(5) 



dnldy^.la,..^.jf.dpl 



are obtained from the interpolated data. The numerical values are given in Table|^and presented 
as functions oixp and y in Fig.[ 
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0.2182 0.58 
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0.3344 1.10 
0.3157 1.60 
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0.57 
0.52 
0.48 
0.38 
0.33 
0.30 
0.31 
0.32 
0.35 
0.37 
0.38 
0.38 
0.38 
0.42 
0.47 
0.62 
0.87 
1.15 
1.60 
2.40 



0.3287 
0.3141 
0.2970 
0.2796 
0.2631 
0.2558 
0.2648 
0.2844 
0.3042 
0.3237 
0.3413 
0.3781 
0.4045 
0.4277 
0.4436 
0.4730 
0.4909 
0.4993 
0.5159 
0.5125 



0.40 
0.31 
0.31 
0.30 
0.30 
0.30 
0.30 
0.30 
0.31 
0.32 
0.32 
0.32 
0.32 
0.33 
0.34 
0.37 
0.43 
0.53 
0.63 
1.00 



0.1543 0.80 
0.1425 0.60 
0.1294 0.55 
0.1168 0.55 
0.1049 0.54 
0.1008 0.54 
0.1062 0.55 
0.1206 0.55 
0.1352 0.58 
0.1506 0.62 
0.1649 0.62 
0.1993 0.62 
0.2268 0.62 
0.2525 0.65 
0.2708 0.70 
0.3053 0.70 
0.3308 0.85 
0.3425 1.10 
0.3664 1.15 
0.3620 2.00 



-0.6 
-0.4 
-0.2 
0.0 
0.2 
0.4 
0.6 
0.8 
1.0 
1.2 
1.4 
1.6 
1.8 
2.0 



0.9795 
0.9994 
1.0120 
1.0010 
0.9796 
0.9373 
0.8941 
0.8448 
0.7898 
0.7291 
0.6611 
0.5807 
0.4942 
0.4069 



0.9311 
0.9521 
0.9688 
0.9646 
0.9317 
0.8751 
0.8159 
0.7597 
0.6905 
0.6131 
0.5296 
0.4451 
0.3612 
0.2806 



Table 6: pr integrated invariant cross section F [mb-c], density distribution dn/dxp, mean 
transverse momentum (px) [GeV/c], mean transverse momentum squared (py) [(GeV/c)^] as 
a function of xp, as well as density distribution dn/dy as a function of y for 7r+ and vr^. The 
statistical uncertainty A for each quantity is given in % 



The Pt integrated -n^ /-n^ ratios, and the first and second moments of the p^ distributions 
are shown as functions of a; in Fig.l27l 

Due to the absence of data with sufficient phase space coverage allowing for integration 
over transverse momentum, a comparison with other experiments is unfortunately not possible. 
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Figure 26: Integrated distributions of vr"*" and n~ produced in p+C interactions at 158 GeV/c: a) 
density distribution dn/dxp as a function of xp', b) invariant cross section F as a function of 
Xi?; c) density distribution dn/dy as a function of y 




[o' 0.8 
> 



■0.6 



0.4 







1 1 1 1 1 
. b) 


1 1 1 1 1 1 1 

•71+ 




. 


:/ 

.1,1, 


.1,1,1, 




-0.1 0.1 0.2 0.3 0.4 0.5 0.6 -0.1 0.1 0.2 0.3 0.4 0.5 0.6 -0.1 0.1 0.2 0.3 0.4 0.5 0.6 



Figure 27: a) 7r+/7r ratio, b) mean px, and c) mean pip as a function of xp for 7r+ and tt 
produced in p-i-C interactions at 158 GeV/c 



8 Dependence on the number of grey protons 

Due to the large fraction of single projectile collisions in minimum bias p-i-C interactions, 
the number of grey protons ngiey measured in this reaction has a steep maximum at zero counts, 
as shown in Fig.|3]above. Correspondingly the event sample decreases rapidly from 377 000 to 
102 000 and 26 000 respectively by selecting events with one or two grey protons. A general 
study of double differential inclusive cross sections as described above is therefore not feasible 
for these subsamples with reasonable statistical errors. On the other hand grey proton selec- 
tion allows an important extension of the physics analysis since it effectively suppresses the 
contribution of peripheral collisions in favour of more central events with multiple projectile 
collisions. 

The problem of limited statistics can be overcome by extracting the pt integrated yields 
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directly from the data sample as a function of xp only. There is a limit to this procedure set by 
the necessities of particle identification. The variation of the total momentum over the chosen 
bin width in a;^ must be small enough to allow for the necessary resolution of the energy loss 
distribution. At low x p this variation becomes prohibitive due to the strong variation of the total 
momentum with pp. The method is therefore limited to xp > 0.025. The quality of the dE/dx 
distributions obtained is shown in Fig.|28lfor two values of xp- 

The extracted pp integrated pion yields are given in Table IT] for the total minimum bias 
sample and for the two centrality selections ngiey > 1 and rigrey > 2. The distributions have been 
corrected by using the treatment developed in Sect. |4| integrated over transverse momentum. 
The yields obtained for the minimum bias event sample are in good statistical agreement with 
the values obtained from the integration of the interpolated double differential cross sections 
given in Tabled 

In comparison to the xp dependence in minimum bias condition, a systematic steep- 
ening is observed as shown in Figs. andl^Db for 7r+ and vr^, respectively. The difference 
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Table 7: dn/dxp for rigrey selection. The uncertainty A is given in % 
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Figure 29: a) dn/dxp for Ug^ey selection and b) the difference 5 between samples with n^^ey 
selection and minimum bias for 7r+. The lines shown are to guide the eye 
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Figure 30: a) dn/dxp for ngrey selection and b) the difference 5 between samples with n^^^y 
selection and minimum bias for tt^. The lines shown are to guide the eye 



5{dn/dxF) = {{dn/dxF)ng,^y — {dn/dxp) mm. bias)/ {dn/dxF)ng,^y is shown separately in Figs.l29b 
andESb. 
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9 Availability of the presented data 

As in lUl the tabulated values of NA49 data are available in numerical form on the Web 
Site |Il6|. In addition two sets of 7r+ and tt^ momentum vectors (5x10'' each) can be found 
on this site. They are generated via Monte Carlo following the data interpolation presented in 
Sect.|5]between the limits -0.1 < xp < 0.5 and < pt <2 GeV/c using a slight extrapolation of 
the data in the backward hemisphere at low px and to pt values beyond 1.8 GeV/c. Normalized 
invariant distributions may be deduced from these vectors using the inelastic cross section given 
in Sect. 13.21 and the following total integrated pion multiplicities in the xp range from -0.1 to 
+0.5: 

(n^+) = 3.279 
(n^-) = 2.909 

10 Conclusions 

A new set of inclusive cross sections on pion production in minimum bias p+C collisions 
at the CERN SPS is presented. The data cover the central production region within the total 
range of -0.1 < Xi? < 0.5 andO < < 1-8 GeV/c for the first time. The statistical uncertainties 
are typically at the few percent level over the 270 measured bins per charge, with systematic 
errors of less than 5%. A detailed discussion of the results, including in particular an in-depth 
comparison to the recently published data from NA49 on p-i-p interactions, will be presented in 
an accompanying paper. 
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